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Abstract 

We eslimated local and metapopulation effective sizes (N e and meta-N e ) for 
three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus 
alpinus) inhabiting a freshwater system comprising seven interconnected lakes. 
First, we hypothesized that N e might be inversely related to within-species 
population divergence as reported in an earlier study (i.e., F ST : S. salar> S. 
fontinalis> S. alpinus). Using the approximate Bayesian computation method 
implemented in ONeSAMP, we found significant differences in N e (NeOSMp) 
between species, consistent with a hierarchy of adult population sizes 
(N e s.s a iar<Nes.fontinaiis<N eS .aipinu S )- Using another method based on a measure 
of linkage disequilibrium (LDNE: NeLDNE ) , we found more finite N e values 
for S. salar than for the other two salmonids, in line with the results above 
that indicate that S. salar exhibits the lowest N e among the three species. 
Considering subpopulations as open to migration (i.e., removing putative 
immigrants) led to only marginal and non-significant changes in N e , suggest- 
ing that migration may be at equilibrium between genetically similar sources. 
Second, we hypothesized that meta-N e might be significantly smaller than the 
sum of local N e s (null model) if gene flow is asymmetric, varies among sub- 
populations, and is driven by common landscape features such as waterfalls. 
One 'bottom-up' or numerical approach that explicitly incorporates variable 
and asymmetric migration rates showed this very pattern, while a number of 
analytical models provided meta-N e estimates that were not significantly dif- 
ferent from the null model or from each other. Our study of three species 
inhabiting a shared environment highlights the importance and utility of dif- 
ferentiating species-specific and landscape effects, not only on dispersal but 
also in the demography of wild populations as assessed through local N e s and 
meta-N e s and their relevance in ecology, evolution and conservation. 
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a Equal effort. 

Introduction 

Multispecies comparisons of closely related taxa in land- 
scape genetics studies have the potential to assist in vali- 
dating models in theoretical population genetics and 
empirically assess their applicability, that is, the extent to 
which they capture the nuances of natural populations 
(Hanfling and Weetman 2006; Fraser et al. 2007; Beebee 



2009). When model results are discordant, conclusions 
can be difficult to draw, but agreement, on the other 
hand, suggests that parameters inferred from the data 
may be robust to differential assumptions characterizing 
the theoretical models (e.g., Whitlock and Barton 1997). 
Studies making these types of comparisons may be espe- 
cially informative for conservation and management of 
wild populations, as they can elucidate the general 
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influence of life history aspects and demography on pop- 
ulation genetics (Turner et al. 1996; Whiteley et al. 2004; 
Manier and Arnold 2006; Gomez-Uchida et al. 2009). 

The effective size of a population (N e , Wright 1931) 
was defined to measure the rate of heterozygosity loss in 
a population, and is explicitly linked to migration rate 
(m) in the classical island model, through F S t (Wright 
1965). In addition to its importance at the interplay of 
population structure and dispersal, N e is a central parameter 
in conservation biology as it is linked to the long-term via- 
bility and extinction risk of populations (Lande 1988; Frank- 
ham 1995). However, only recently have methods been 
developed that use genetic data to assess N e independently of 
m (reviewed Wang 2005; Luikart et al. 2010). Such methods 
have been applied in numerous studies, mostly to reveal pat- 
terns of variable N e within many taxa (fishes: Fraser et al. 
2004; Gomez-Uchida et al. 2008; Palstra et al. 2009; rep- 
tiles: Manier and Arnold 2005; mammals: Kaeuffer et al. 
2007; birds: Hoeck et al. 2010; see Palstra and Ruzzante 
2008 for review). On the other hand, comparisons of N e 
among different species occupying the same habitat have 
been sparse, even though these studies may provide evo- 
lutionary insights into previously documented ecological 
interactions (e.g., competition: Manier and Arnold 2005; 
predator-prey: Manier and Arnold 2006). Such analyses 
may also elucidate the role of life history in shaping vari- 
ation in Ne in a shared common environment. 

Most species naturally occur in networks of more or 
less interconnected populations (hereafter called meta- 
populations), thereby potentially biasing methods for 
estimating N e that assume single closed populations. 
More importantly, this currently limits our understanding 
of the magnitude and behavior of N e under typical condi- 
tions in the wild, both at the local scale of a population 
and at the global scale of a meta-population. Solutions 
have been proposed at different levels with various 
approaches (e.g., Fraser et al. 2007; Palstra and Ruzzante 
2008; Waples and Do 2010). At the population level, 
approaches that simultaneously calculate contemporary 
N e and m have gained popularity (Vitalis and Couvet 
2001; Wang and Whitlock 2003; but see Waples and Eng- 
land 2011). Alternatively, some authors have attempted to 
identify immigrants and assess how they affect temporal 
estimates of contemporary N e (Walter et al. 2009). Very 
few studies so far have explored the effects of gene flow 
on single-sample estimators of contemporary N e (e.g., 
Palstra and Ruzzante 2011; Waples and England 2011). At 
the metapopulation level, the concept of effective meta- 
population size (meta-Ne), that is the effective size of a 
subdivided population, has been the focus of extensive 
theoretical work (e.g., Whitlock and Barton 1997; Nunney 
1999). Wang and Caballero (1999) argued that meta-Ne 
can be smaller or larger than the sum of N e , depending 



on whether all subpopulations have equal variance in 
reproductive success. Recent empirical studies on salmo- 
nid fish have reported reductions in meta-N e in the 
context of asymmetric gene flow and variable local N e 
(Kuparinen et al. 2010; Palstra and Ruzzante 2011). These 
findings hence suggest that meta-N e may indeed often be 
much smaller than would be predicted in relatively simple 
metapopulation models, but clearly more empirical esti- 
mates are required to validate this generalization. 

This study focuses on three salmonid fish species, 
Atlantic salmon (Salmo salar), brook trout (Salvelinus 
fontinalis), and Arctic charr (Salvelinus alpinus), living in 
sympatry in a system of seven interconnected lakes in the 
Upper Humber River, Gros Morne National Park (New- 
foundland, Canada: Fig. 1). All three species form discrete 
spawning aggregations within each lake, but with varying 
degrees of global population divergence (P ST : S. salar > S. 
fontinalis > S. alpinus; Table 1; Gomez-Uchida et al. 
2009). Such levels of population divergence seem to be 
inversely related to population density as local mark- 
recapture studies indicate the abundance of S. alpinus > 
S. fontinalis > S. salar (Table 1; see also Ryan and Kerekes 
1988; Anions 1994). 

Our goal in the present paper is twofold. We first test 
the hypothesis that contemporary effective population size 
will be inversely related to the degree of population struc- 
ture with S. alpinus and S. salar exhibiting, respectively, 
the highest and lowest estimates of effective population 
size (N e ) with S. fontinalis exhibiting an intermediate N e . 
Secondly, we estimate meta-Ne (meta-N e ) using a variety 
of models and then compare these estimates across mod- 
els within species. Our second hypothesis, now related to 
meta-Ne, is that meta-N e can be significantly smaller than 
the sum of local N e s if gene flow is asymmetric and varies 
among subpopulations (Tufto and Hindar 2003). 

Materials and Methods 

Study system and genetic data 

The focal system consists of seven interconnected lakes or 
'ponds' (P1-P7: Fig. 1) in the Upper Humber River, Gros 
Morne National Park, Newfoundland (details in Gomez- 
Uchida et al. 2009), where Atlantic salmon (Salmo salar), 
brook trout (Salvelinus fontinalis) and Arctic char (Salv- 
elinus alpinus) naturally co-occur as non-anadromous 
landlocked populations, probably founded from anadro- 
mous populations that survived the Wisconsinan glacial 
maximum (10,000-12,000 BP) (Batterson and Catto 2001; 
Shaw et al. 2006). 

The genetic data correspond to a subset of these ponds 
(hereby subpopulations) (Gomez-Uchida et al. 2009) 
including individual samples from S. salar, S. fontinalis 
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Figure 1. (a) Study site in the Upper Humber River, Gros Morne National Park (rectangle shows the area in western Newfoundland, Canada), (b) 
The population system: seven interconnected lakes locally referred to as "ponds" (P1-P7); P1 and P7 are found above waterfalls (in blue) of 5 m 
and 3 m, respectively, which are impassable for fish, whereas P2-P6 are found below waterfalls. River flows naturally from P1 to P5 and from P7 
to P5. (c) Estimates of migration rate for (a) Salmo salar, (b) Salvelinus fontinalis, and (c) Salvelinus alpinus from a previous study (Gomez-Uchida 
et al. 2009) showed that these waterfalls represent barriers to dispersal and probably cause asymmetrical gene flow (arrows in log-scale), even 
though the magnitude was species-specific. Red lines indicate that no estimates of gene flow could be reliably obtained because samples were 
unavailable or insufficient (see text for details). 



Table 1. Data set and population (genetic and ecological) parameters for a community of salmonids {Salmo salar, Salvelinus fontinalis, and Salv- 
elinus alpinus) from the Upper Humber River, Gros Morne National Park (Newfoundland, Canada). 



S. salar 



S. fontinalis 



S. alpinus 



Remarks 



Data set 
Number of genotyped 

microsatellite loci 
Mean number of alleles 

per locus 
Number of subpopulations 

or ponds (names) 
Number of genotyped 

individuals 
Population 
Global population 

divergence - F' ST (95% CI) 
Global inbreeding 

coefficient - F' n 
Average immigration 

rate (m± SE) 
Variance of dispersion 

distance (<r 2 ) 
Adult population size 



12 
13.2 

6 (P1-P6) 
456 

0.202 (0.200-0.205) 
0.202 

0.053 ± 0.025 

0.026 

989 



13 
8.6 

7 (P1-P7) 
706 



11 

18.2 

3 (P2-P4) 
301 



Details in Table S1-S4 from 

Gomez-Uchida et al. (2009) 
Over all loci and populations 

Minimum sample size of 20 to 

estimate N e 
Maximum of three loci with missing 

genotypes per individual* 



0.076 (0.058-0.089) 0.020 (0.012-0.025) Standardized by heterozygosity 



0.088 

0.061 ± 0.034 

0.023 

2021 



0.012 

0.153 ± 0.050 



12773 



Standardized by heterozygosity 

Contemporary estimator using BIMr 

No isolation by distance in S. alpinus 

Mark-recapture estimate from 
Hardings pond (P2) 



♦Requirement for ONeSAMP (Tallmon et al. 2008). 



and S. alpinus genotyped for 12, 13, and 11 microsatellite 
loci, respectively (Table 1). Most microsatellites assayed 
for S. salar and S. fontinalis were species-specific, but the 
majority of microsatellites assayed for S. alpinus originated 



from related salmonids (see Supporting Information in 
Gomez-Uchida et al. 2009). However, a high number of 
alleles per locus for S. alpinus, compared with the other 
two species (Table 1), suggest that ascertainment bias is 
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probably negligible, if present at all. We set a minimum 
size of 20 individuals per population sample to estimate 
N e and meta-N e , resulting in the exclusion of one 
subpopulation (P7) for S. salar, none for S. fontinalis, and 
two (P5 and P7) for S. alpinus and leaving 6, 7 and 3 
(total 16) subpopulations available for analyses for each 
species, respectively (Table 1). Population genetics param- 
eter estimates required for some models of meta-N e , 
namely global estimates of divergence and inbreeding and 
contemporary estimates of migration, were taken from 
Gomez-Uchida et al. (2009). 

Effective population size (JV e ) 

Firstly, contemporary N e was estimated for each local 
subpopulation and species using ONeSAMP 1.2 (hereafter 
N e o SM p ; Tallmon et al. 2008). Briefly, the program 
employs approximate Bayesian computation to estimate 
the size of an ideal Wright-Fisher population (with no 
migration or selection) based on summary statistics calcu- 
lated from the empirical data. The timescale to which the 
resulting N e estimate applies is not explicit, but probably 
applies to recent generations (Luikart et al. 2010). We 
obtained four replicate estimates of N e and their 95% 
credible limits (CL) under the following prior N e values: 
NeMiN = 10 and N eMAX = 2000. Although N eOSM p 
appeared robust to different N eMW values across subpop- 
ulations and species (e.g., N eM i N = 20, 40: authors' un- 
publ. results), changes in N eMAX were not assessed, also 
because effective sizes in the three salmonids species were 
initially thought unlikely to be much larger than N eMAX 
given the size of these small ponds. 

Secondly, we estimated N e plus parametric 95% confi- 
dence intervals (CI) using the linkage disequilibrium 
(LD) of LDNE (hereafter N eLDNE . Waples and Do 2008). 
LDNE estimates reflect the number of parents that con- 
tributed to the sample (Waples 2005) and assume that 
LD at unlinked loci arises solely from genetic drift in an 
isolated population (Hill 1981; Wang 2005). Following 
Waples and Do (2010), we set the cut-off probability 
(P crit ) for alleles of low frequency to be included in the 
estimation at P crit = 1/2S, where S is the subpopulation 
sample size. 

Effects of immigration on A/ eLDNE and N eOSmP 
assuming open subpopulations 

The estimator methods implemented here assume that 
populations are closed to migration (Tallmon et al. 2008; 
Waples and Do 2008), but violation of this assumption 
may result in estimation bias (Vitalis and Couvet 2001; 
Palstra and Ruzzante 2011; Waples and England 2011). 
We thus conducted exploratory analyses to assess the 



impact of immigration on N e , by identifying and exclud- 
ing putative migrants from samples, and repeating the 
estimation of N e in ONeSAMP and LDNE. 

First-generation immigrants were identified in GENE- 
CLASS 2 (Piry et al. 2004) through Bayesian assignment 
of genotype likelihoods (Rannala and Mountain 1997). 
Here, individual genotypes were ranked according to the 
ratio Lh/imax that relates the likelihood of drawing geno- 
types from the populations in which they were sampled 
(numerator) with the maximum likelihood of such geno- 
types considering any of the study populations (denomi- 
nator). We assumed that all sources of migrants were 
sampled (Paetkau et al. 2004) and attempted to minimize 
erroneous identification of migrants by setting the type I 
error at a = 0.01 and by graphically assessing likelihood 
scatter plots between recipient and source populations. 

Effective metapopulation size (meta-/V e ) 

The estimation of meta-N e from genetic data depends on 
the underlying model of spatial genetic structure (Wang 
and Caballero 1999). We have assumed that within our 
system, extinction and recolonization were infrequent and 
thus our model was more related to a "patchy" or "main- 
land-island" metapopulation than a classic Levins meta- 
population (Koizumi et al. 2006). Here, we estimated 
meta-N e assuming seven different models (Table 2), 
Firstly, we simply summed the various estimates of local 
N e (E(AT e )); this is our Null model. We then estimated the 
associated 95% confidence intervals following the standard 
equation SE(E (N e ) = sqrt [E (SE(N e ) 2 )]. This was done 
taking into consideration the fact that lower and upper 
confidence intervals differed from each other. Then, we 
estimated meta-N e as defined by the Island (Wright 1931), 
Stepping Stone (Maruyama 1970), and Neighborhood 
models (Wright 1946; Maruyama 1971). These models 
assume that all subpopulations contribute equally to the 
next generation, but differ regarding assumptions on the 
spatial configuration and connectivity of subpopulations, 
yielding meta-N e that should be larger than Null estimates 
(Table 2). Population parameters included in these models 
were taken from Table 1 or indirectly calculated from data 
obtained in Gomez-Uchida et al. (2009). Effective popula- 
tion density (D) for Neighborhood meta-N e was obtained 
by dividing the sum of subpopulation N e s by the linear 
habitat length (L = 11 km); the variance of dispersion dis- 
tance (cr 2 ) was estimated from the slope (b) of the isola- 
tion-by-distance relationship described in (Gomez-Uchida 
et al. 2009) according to 4Der 2 = lib (Rousset 1997). Step- 
ping Stone and Neighborhood meta-Ne were ignored for 
S. alpinus, because we found no indications that an 
isolation-by-distance model was applicable to this species 
(Gomez-Uchida et al. 2009). Second, we calculated 
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Table 2. Models of spatial genetic structure and estimation of effective metapopulation size (meta-W e ) or the size of an idealized population with 
the same rate of inbreeding observed in the subdivided population under study. 



Name 



Model* 



Expectations 



Reference 



Island 


meta 


-We 


/=! 
~~ 1-% 




s 

meta-We > EWe, unless ^57= 0; inequality 

/=1 

IMLI CoSCj W I LI I 15J ^UIVtrlUtrMLtry 


Wright (1943) 


Stepping stone 
(circular) 


meta 


-W e 


= E« ei 

1=1 


+ ^=_ 


5 

meta-We > £We,; inequality increases with 

/=1 

decreasing average immigration rate (m) 


Maruyama (1970) 


Neighborhood 
(linear) 


meta 


-W e 


= E«e, 
1=1 




5 

INfcrLd /V e / . /Ve; 1 IMeqUdllLy IMLIcdbfcrS VVILM 
1=1 

increasing length of the habitat (A.), but 
decreases with increasing population density (D) 


Wright (1946) 
Maruyama (1971) 


Interdemic 


meta 


-We 






meta-We < Z^e,; inequality increases with 

;=1 

metapopulation inbreeding 


Nunney (1999) 


Spatiotemporal 


meta 


-We 


= E«e, 
1=1 


(i+w) 


meta-We > Z^; inequality increases with 

1=1 

the average number of migrants per 
generation (M = W e m) 


Kobayashi and Yamamura 
(2007) 



*Metapopulation parameters: s = number of subpopulations/ponds; W e = average subpopulation W e ; F ST = global genetic divergence among sub- 
populations; m= average immigration rate; D = linear population density (individuals/km); L = length of habitat (km); a 2 = variance of dispersion 
distance; F| T = coefficient of global (metapopulation) inbreeding. 



Interdemic meta-N e ; this model predicts meta-N e can be 
smaller than the Null estimate when accounting for vari- 
able fitness among subpopulations (Nunney 1999). Third, 
we explored Spatio-temporal meta-N e by Kobayashi and 
Yamamura (2007), which provides similar expectations to 
Island meta-N e , but is centered on the average number of 
contemporary immigrants (Table 2). Finally, we used the 
empirical 'bottom-up' approach of Tufto and Hindar 
(2003), which combines subpopulation N e s (Ne^calcu- 
lated previously using ONeSAMP) with unidirectional 
migration rates (estimated using BIMr) to calculate T&H 
meta-N,,. Here, we minimized the eigen-value of the 
resulting metapopulation matrix employing an R package 
(R Development Core Team 2010) developed by J. Tufto 
(available from http://www.math.ntnu.no/~jarlet/migra- 
tion/). Given sample size limitations (see above), meta-N e 
for S. fontinalis was estimated considering all seven 
subpopulations and for S.salar, it was estimated over six 
subpopulations (omitting P7). Salvelinus alpinus was 
collected in adequate numbers only in three ponds (P2-P4); 
therefore, meta-N e for this species was limited to those 
three ponds (see Results). 

Statistical analyses and hypothesis testing 

We assessed the null hypothesis of no correlation between 
N e and sample size using Kendall's t correlation coeffi- 
cient and two-tailed exact probabilities, given the uncer- 
tainty as to whether the data were drawn from a bivariate 
normal distribution or not. Such correlations permitted 



to evaluate any dependence of N e on sample size using 
ONeSAMP as LDNE provides unbiased estimates (Waples 
and Do 2008). To evaluate differences in N e estimates 
among species, the null hypothesis of no differences was 
tested using nonparametric Kruskall-Wallis tests for mul- 
tiple species comparisons and Mann-Whitney one-tailed 
tests for pairwise species comparisons. When replicates 
were available (e.g., ONeSAMP), we compared the har- 
monic means to minimize the contribution of overly dis- 
tinct replicates. All tests were implemented using R 2.13.2 
(R Development Core Team 2010). 

Results 

W e : Comparisons between estimators and 
species 

N e was estimated for 16 species/subpopulation combina- 
tions (Fig. 1). While all 16 N e osMP values were finite, their 
magnitude differed across replicates, in some cases con- 
siderably so (Fig. 2a-c). For N e LDNE, on the other hand, 
only seven or eight estimates were finite depending on 
whether immigrants were included or excluded, most of 
which were concentrated among S. salar subpopulations 
(Fig. 2d-f). The range of values (box plots) for N c osmp 
replicates was narrow for S. salar, intermediate for S. fon- 
tinalis and wide for S. alpinus (Fig. 2a-c). Two N c osmp 
replicates exhibited extremely high upper confidence lim- 
its (in the millions of individuals). On the other hand, 
only four upper 95% CIs were finite with LDNE (Fig. 2). 
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S. salar 



S. fontinalis 



S. alpinus 




(b) 





(e) 



P1 P2 P3 P4 P5 P6 P1 P2 P3 P4 P5 P6 P7 P2 P3 P4 



Subpopulation 



Figure 2. Estimates of effective population size (W e ) and their uncertainty for Salmo salar (left column), Salvelinus fontinalis (middle column), and 
Salvelinus alpinus (right column) using two estimators, N,,osmp (top row; boxplots) and A/ etDNf (second row; empty circles). Blue boxplots or filled 
circles are estimates of A/ c following exclusion of immigrant genotypes (see text for details), thus assuming populations open to gene flow. Red 
dots in boxplots are the harmonic means among four N^osmp replicates, whereas green triangles represent the lowest (inverted triangles) and 
highest 95% CIs among four N e osMP replicates. Missing or out-of-bound empty or filled symbols, upper 95% CIs, or both, suggest infinitely large 
values. 



Finite N e estimates were generally poorly correlated with 
sample size (S. salar, N e osMP t = 0.69, P = 0.06; S. salar, 
Ncldne T = -0.10, P = 0.80; S. fontinalis, N e osMP 
t = —0.14, P = 0.77; no correlation test was performed 
with N e LDNE for S. fontinalis and no correlation tests were 
attempted for S. alpinus with only three estimates avail- 
able). Comparisons of N e osMP harmonic means among 
the three salmonids indicated significant differences 
(Kruskal-Wallis yf = 8.8, df = 2, P = 0.01); furthermore, 
one-tailed Mann-Whitney pairwise tests were consistent 
with a hierarchy of N e osMP values as proposed (S. salar < 
S. fontinalis: P < 0.05; S. fontinalis < S. alpinus: 
P < 0.008; S. salar < S. alpinus: P < 0.01). Complete 
details on sample sizes, N e iDNE> and N e osMP can be found 
in Table SI (Supporting Information). 

Immigrant genotypes and their impact on 
Weosivip and N eLD NE (open subpopulations) 

Immigrant genotypes that could be identified ranged 
between 1 and 3 among subpopulations and across spe- 
cies (Table SI). In four out of five S. fontinalis 
subpopulations, the likely source of immigrants was 
either PI or P7, both ponds found above waterfalls 
(Fig. 1), whereas all three S. salar individuals identified as 



immigrants in P3 probably originated from P5 (Table 
SI). The weak population structure observed in S. alpinus 
prevented the reliable identification of potential immi- 
grants for this species. 

Both estimators - N e0 sMP and NeLDNE - suggested that 
the changes in N e after exclusion of immigrants were neg- 
ligible and non-significant (Fig. 2, Table SI). 

Meta- N e 

Given the generally low precision of NeLDNE. the 
estimates of meta-N e were based only on NeosMP. Meta- 
Ne differed across species by up to three orders of 
magnitude; it was smallest in S. salar, intermediate in S. 
fontinalis, and largest in S. alpinus. Moreover, these quali- 
tative differences across species were constant regardless 
of the spatial population model assumed (Fig. 3). Indeed, 
the differences in meta-N e between estimator models were 
very small, with one notable exception: the meta-N e 
obtained with the Tufto and Hindar (2003) model, which 
were always much smaller than any of the other meta-N e 
estimates. In fact, T&H meta-N e for S. salar and 
S. fontinalis were similar to the N e0 sMP estimates 
obtained for both species in one subpopulation, PI (Table 
SI and Table 3). 
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S. salar 




1 2 3 4 5 6 7 



S. fontinalis 



(b) 



1 2 3 4 5 6 7 

Model No. 



S. alpinus 



• Null 
o Island 

□ Stepping-stone 
o Neighborhood 
A Interdemic 
V Spatiotemporal 
x Tufto-Hindar 



2 3 



(c) 



5 6 7 



Figure 3. Estimates of mefa- N c from the null model (sum of subpopulation A/ e 's) and six metapopulation models (island, stepping stone, 
neighborhood, interdemic, spatiotemporal, and Tufto-Hindar: see legend) for (a) Salmo salar, (b) Salvelinus fontinalis, and (c) Salvelinus alpinus. 
Stepping stone and neighborhood models require the assumption of isolation-by-distance and were therefore omitted for S. alpinus (see text for 
details). 



Discussion 

In this study, we estimated the local and metapopulation 
effective sizes of three salmonid fish species living in 
sympatry in a spatially fragmented system of seven inter- 
connected subpopulations and exhibiting marked differ- 
ences in demographic attributes as measured by the F sx 
analysis of Gomez-Uchida et al. (2009). First, putative 
immigrants had negligible effects on estimates of local 
effective population size (N e ) and the estimates generally 
fit our expectations based on demography and life history 
with the local effective size of S. salar < S. fontinalis < S. 
alpinus. Second, we implemented six approaches to esti- 
mate effective metapopulation size (meta-_iV e ) and com- 
pared these estimates with the sum of local N e J^(N e ), a 
model we defined as Null. The majority of meta-Ne esti- 
mators yielded results indistinguishable from the sum of 
local sizes, despite considerable variation in population 
structure among the species. This suggests that these ana- 
lytical models may not fully capture the complexity of 
natural systems or simply are not well suited for applica- 
tion to empirical data. On the other hand, the estimates 
obtained with the Tufto and Hindar (2003) approach 
[T&H meta-Ne] differed substantially from the other mod- 
els and were strongly influenced by connectivity patterns 
inferred from the genetic data. The linear property of our 
system, combined with asymmetric dispersal facilitated by 
the presence of waterfalls, may thus have resulted in strong 
reductions in meta-N e in line with other recent empirical 
work and theoretical expectations. Below, we review some 
caveats regarding the methods used to estimate contempo- 
rary local N e and discuss how immigration seemed to have 
negligible impacts on N e ; we then argue how possible dif- 



ferences in life history (migration and abundance) reflect 
on differences in average local N e s among species. We sub- 
sequently consider estimates of meta-Ne and the similari- 
ties and differences among models within species. We 
finish by explicitly discussing our results in the context of 
known species differences in life history (migration as 
measured by gene flow) and abundance. 

ONeSAMP and LDNE: Caveats and patterns 
across species 

The two methods differed in performance within the scope 
of our study. Firstly, they differed in the number of finite 
N e estimates they provided: all point estimates of effective 
population size were finite when estimated with ONeS- 
AMP, but only seven or eight (depending on whether 
immigrants were included or excluded) were finite when 
estimated with LDNE. As the estimation of contemporary 
N e depends largely on the genetic drift signal, which is 
stronger in smaller populations (Fraser et al. 2007; Palstra 
and Ruzzante 2008; Waples and Do 2010), it should not be 
surprising that both approaches provided finite values only 
for S. salar, the species with the lowest abundance and his- 
torical census size (Ryan and Kerekes 1988; Anions 1994). 

Secondly, the two methods apply to different time 
scales: while LDNE uses linkage disequilibrium to esti- 
mate the effective size in the parental generation (Waples 
and Do 2008), ONeSAMP uses several metrics, including 
linkage disequilibrium (Tallmon et al. 2008) to estimate 
effective size over the last few, but undefined number of 
generations (Luikart et al. 2010; see also Skrbinsek et al. 
2012). This is the likely cause of the discrepancy of results 
for S. salar in the headwater pond PI (Fig. 2 and Table 
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SI): the M-ratio test (Garza and Williamson (2001), 
which is considered in ONeSAMP, suggests that a bottle- 
neck has occurred in this population (Gomez-Uchida 
et al. 2009), but this test statistic is not considered in the 
LDNE method. It should be noted, however, that the 
removal of first-generation migrants prior to the estima- 
tion of effective size does not necessarily completely elim- 
inate the influence of migration if migration took place 
in previous generations. Such earlier migration could be 
affecting N e particularly when estimated with ONeSAMP. 

A third caveat regarding our replicated estimates of N e 
is that they were all obtained with priors set at 
N em in = 10 and N emax = 2000. For S. salar, only few point 
estimates (3 out of 24) exceeded 2000 (Table SI), suggest- 
ing that the true N e s for this species falls mostly within 
this range. On the other hand, approximately one-third 
of the replicates for S. fontinalis, and two-thirds of the 
replicates for S. alpinus were higher than the N max prior 
(See Table SI). These results suggest that for these two spe- 
cies, and particularly for S. alpinus, the true N e values may 
indeed be larger than 2000 and hence the program had dif- 
ficulty estimating them (David Tallmon, personal commu- 
nication). Our caveats notwithstanding, Beebee (2009) also 
found that N e0 sMP showed higher precision than N e LDNE 
for natterjack toad (Bufo calamita) populations for which 
N varied between 2 and 500, that is, within the range of the 
estimates found for S. salar in our study. 

Additionally, age-structure and overlapping generations 
can significantly complicate the estimation of N e (Palstra 
et al. 2009). We used samples consisting of different ages, 
which probably reflect a quantity intermediate between 
N e and Nb, the annual number of breeders (Waples 
2006). The samples analyzed in this study may include 
three or four cohorts based on length ranges (S. salar. 78 
-248 mm; S. fontinalis: 50-230 mm; S. alpinus: 87- 
157 mm, authors' unpublished results) and length-at-age 
data from earlier surveys in Gros Morne National Park 
(Rombough et al. 1978; Barbour et al. 1979; McCairns 
2004). The number of cohorts in our samples thus 
approximates the mean generation length (G) based on 
estimates of the mean age of reproductive parents in this 
system (S. salar G = 4.5: (Rombough et al. 1978); S. fon- 
tinalis G = 4.5: (McCairns 2004); S. alpinus G = 3.5: 
(Barbour et al. 1979). Under such sampling conditions, 
Waples and Do (2010) suggested that Nb « N e , even 
though other researchers have opted to report their esti- 
mates as Nb (e.g., (Beebee 2009) or apply the correction 
N e = GNb (e.g.) Ardren and Kapuscinski 2003; see also 
Robinson and Moyer 2012). Therefore, it may be recom- 
mendable to complement single- with two-sample tempo- 
ral estimates, especially methods tailored to gauge allele 
frequency shifts between cohorts, which require individual 
ages (Palstra et al. 2009; Waples and Do 2010). 



To summarize, the local effective sizes for the S. alpinus 
and S. fontinalis populations studied here are likely larger 
than can reliably be estimated given the sample sizes (S~ 
90-100) and numbers of microsatellite loci (L-12-13) 
available, but this limitation appears to be stronger for 
LDNE than for ONeSAMP. Overall, N e0 sMP may there- 
fore appear to be more general than N eLDNE , due to the 
rationale and flexibility of its methodology (prior specifi- 
cations, simulation approach). On the other hand, the 
main drawback of ONeSAMP may be its undefined tem- 
poral reference (Luikart et al. 2010; see also Skrbinsek 
et al. 2012), thereby reducing its utility for genetic moni- 
toring. Ultimate choice of N e estimator should be guided 
by the specific aims and goals of each study. 

Impacts on /V e osMP and A/ eLDNE of accounting 
for gene flow 

We attempted to account for the effect of gene flow 
among subpopulations by re-estimating effective popula- 
tion sizes after eliminating putative immigrants. This 
approach was implemented with the two species for 
which immigrants were identified: S. salar and S. fontinal- 
is. Although in all cases the point estimates of N e 
decreased after the removal of putative migrants, these 
decreases were minor and non-significant (Fig. 2), sug- 
gesting that migration may be sustained and at equilib- 
rium between genetically similar sources (Waples and 
England 2011). We reiterate, however, that the removal of 
first-generation migrants does not necessarily completely 
remove the influence of migrants in the estimation of N e . 
As linkage disequilibrium between physically unlinked loci 
decreases by 50% each generation, any LD caused by 
migration in the grand-parental or great-grand-parental 
generation would probably still influence the estimation 
of N e with either method. 

Comparisons of N e among species 

Despite the limitations imposed by the relatively large 
population sizes for two of our study species, both 
approaches provided similar insight into putative links 
between population abundance, life history and N e . Pair- 
wise comparisons highlighted that S. alpinus probably has 
higher N e than the other two species, matching expecta- 
tions based on its life history and abundance (e.g., high 
population abundance and connectivity (Gomez-Uchida 
et al. 2009). Landlocked S. alpinus from Gros Morne 
National Park aggregate in large numbers in only three 
subpopulations of the watershed (Gomez-Uchida et al. 
2009) and are rare or even absent in other lakes within 
the system (authors' personal observations). Arctic charr 
were described as having slow growth, low fecundity, and 
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early age at maturity by previous authors (3 years, Romb- 
ough et al. 1978), all likely consequences of intense inter- 
and intraspecific competition in low-productive acidic 
lakes (Rombough et al. 1978; Anions 1994). The adult 
census estimate of N ~ 13,000 for S. alpinus within one 
pond (P2) based on mark-recapture (Anions 1994) sug- 
gests that N e for this species should indeed be large, at 
least in the order of a few thousands, even though impre- 
cision of this effective population size estimate may be an 
issue. Massive densities of so-called dwarfed or 'stunted' 
S. alpinus populations have also been described in Norwe- 
gian oligotrophic lakes (Jansen et al. 2002; Persson et al. 
2007). Conversely, the estimated adult census size for S. 
salar within P2 was N ~ 1000 (Anions 1994), a number 
that is consistent with the estimate of effective size 
obtained for this species within this pond (Fig. 2). We 
found significant differences in Nmsmp between S. salar 
and S. fontinalis, two salmonids that showed a similar 
range of distribution in the watershed (in all subpopula- 
tions but one, P7) and comparable average immigration 
rates, but a nearly twofold difference in adult population 
sizes, a finding that is consistent with our hypothesis. 

Effective metapopulation size (meta-JV e ) 

Theory predicts that, under the simplified and restrictive 
assumptions of the island model, the effective size of a 
metapopulation will be higher than the sum of subpopula- 
tion Ne (reviewed in Wang and Caballero 1999). Unsur- 
prisingly, these predictions are borne out in the analyses 
that are explicitly underpinned by those assumptions. The 
difference between (N e ) and island model meta-N e was 
proportional to the degree of fragmentation, with S. salar 
exhibiting the most pronounced spatial structure (highest 
F ST ), and also the largest difference between (N e ) and 
meta-Ne. On the other hand, we find a remarkable consis- 
tency across estimates of meta-N,. for the majority of ana- 
lytical models considered, despite the considerable range of 
population structuring observed (0.020 < F ST < 0.202). 
Moreover, within species, these models provided highly 
similar estimates of meta-N e that did not deviate strongly 
from ^(N e ). Although these models make very different 
assumptions about the metapopulation and consider dif- 
ferent genetic parameters, in practice, their meta-N e esti- 
mates can often be reasonably approximated by simply 
summing individual population N e estimates. Further test- 
ing of these models in different systems and species is nat- 
urally required to assess the validity of this observation. 

The notable exception to this pattern was the meta-Ne 
estimate obtained with the empirical method of Tufto 
and Hindar (2003). This approach provided estimates, 
which were considerably smaller than ^(N e ) as expected 
under conditions of strongly asymmetric gene flow 
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among metapopulation components (Whitlock and Bar- 
ton 1997; Nunney 1999). Indeed, the Tufto and Hindar 
(2003) model makes fewer restrictive assumptions about 
metapopulation dynamics than the other methods, allow- 
ing for consideration of heterogeneity in the size of meta- 
population components and gene flow among these 
components. For instance, in the present study, both S. 
salar and S. fontinalis were very abundant in one of the 
headwater ponds or subpopulations (PI), which lie above 
a waterfall that impedes upstream migration. This is also 
the case for P7, another headwater subpopulation where 
S. fontinalis was found in great numbers. Under such 
conditions, the reproductive contributions among demes 
may diverge strongly from expectations of equality. 
Indeed, for both species, meta-N e was very similar to the 
estimate of local population effective size in the upstream 
(headwater) population. Hence, genetic diversity in this 
(linear) system, for at least two study species, may be lar- 
gely controlled by the diversity that exists in the 
'upstream' source population. This inference is compati- 
ble with Pringle et al. (2011), who suggested that in sys- 
tems with asymmetric gene flow, the upstream edge can 
influence genetic diversity throughout the entire species' 
range. Conversely, in S. alpinus, we documented popula- 
tions with typically high abundance (Table SI) and rela- 
tively high and symmetric migration patterns among 
populations, and meta-Ne was intermediate between the 
smallest and largest subpopulations N e . Therefore, our 
results provide further empirical support for the theoreti- 
cal expectation (Whitlock and Barton 1997) that strong 
inequality among demes in contribution to the shared 
gene pool can reduce the effective size of a subdivided 
population. 

To synthesize, in this study, we have estimated local or 
subpopulation effective sizes, as well as gene flow among 
them, in three sympatric species of salmonids fishes. We 
then compared the sum of these local effective sizes 
^(N e ) to the meta-Ne obtained under the Tufto and 
Hindar (2003) model. The relative magnitude of T&H 
meta-Ne across species was consistent with estimates of 
relative abundance available for one of the ponds at a dif- 
ferent point in time (Ryan and Kerekes 1988; Anions 
1994). Moreover, this pattern mirrors the hierarchy and 
degree of asymmetry in population connectivity described 
in Gomez-Uchida et al. (2009). As our study species 
shared a common environment, differential riverscape 
effects are unlikely to have shaped these patterns. Rather, 
the apparent relation between connectivity patterns and 
census population sizes across species which share a com- 
mon environment, provides further support for the 
notion that density- dependent dispersal may play an 
important role in the maintenance of genetic diversity in 
aquatic systems and in metapopulations in general. 
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